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O . Abstract. We test a set of lattice gauge actions for QCD that suppress small plaquette values and 

^ ' in this way also suppress transitions between topological sectors. This is well suited for simulations 

. in the e-regime and it is expected to help in numerical simulations with dynamical quarks. 

(N : 
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I Our aim is to study the possibility of simulating lattice QCD with a gauge action that 

^ ■ strongly reduces the occurrence of small plaquette values. A gauge background with 

! such a feature is expected to improve the locality properties Jll] of the Overlap-Dirac 

\ operator Dqv 12] . By the same argument one also expects to ease the numerical evaluation 

^ ■ of Dov itself, and help in general dynamical simulations. It can be proven JlLlsl] that as 

I ! long as all plaquette values in a gauge configuration satisfy: 5p := 1 — ^Re Tx{Up) < 

^ I 1 /20.5 then no change of the topological sector is possible. Hence a suppression of low 

■ plaquette values entails a suppression of 2top changes. 

. Simulations constrained in a fixed topological sector can be problematic for evaluating 

physical observables in QCD, where all sectors have to be taken into account with 
^ the correct weight. However such a constraint is perfectly suited for studying QCD in 

^ . the e-regime where predictions exist for observables defined in fixed topological 

j> ! sectors o], which turns the limitation mentioned above into an advantage. However 

ITji \ there are some conditions. The physical volume should be at least L > 1.1 fm |0,|7l]. 

- ■ Moreover, in order to reach very small pion masses, the chiral properties of the Dirac 

operator are crucial. Finally a sound definition of 2top is important to compare with 
predictions in fixed topological sectors. All of these requirements are provided by 
Ginsparg-Wilson fermions. They have an exact, lattice modified chiral symmetry fsl], 
and the fermionic index defines Qtop ^ - Results in the e— regime with Ginsparg-Wilson 
fermions were obtained for the Dirac spectrum ||^ lUJ and for meson correlation 



function s iFtL [12L [TbL [3, which were compared with quenched Chiral Perturbation 
Theory jHMlnr 

Simple examples of gauge actions that suppress small plaquette values (still expected 
to be in the same universality class as Sp) are 

liSf„{Up) = /3- if 5p<e, and +00 otherwise (1) 

' [1 — £ op) 

PST,n{Up) = l5Sp + e-'S'^ (2) 
pS7JUp) = l5Spoxp[e-'s"p] (3) 



TABLE 1. Results for the S'2'\ on a 16^ lattice for various e ' and j3. Pw Is the coupling resulting in 
the same ro with the Wilson action. 



£ ' i3 ro/a 






rP'^'^iM fj 


dt 


Acceptance 



1 


6.18 1 7.14(3) 1 6.18 


7(1) 


7(1) 1 0.015 1 0.1 1 >99% 1 


1 1.00 


1.5 1 6.6(2) 1 6.13(2) 


2.0(1) 


n.a. 1 0.0027 | 0.05 | > 99% | 


1 1.18 


1.0 1 7.2(2) 1 6.18(2) 


1.3(1) 


7(1) 1 0.0014 1 0.02 - 0.01 1 >99% | 


1 1.25 


0.8 1 7.0(1) 1 6.17(1) 


1.1(1) 


9(1) 1 0.0025 1 0.1 1 >99% | 


1 1.52 


0.3 1 7.3(4) 1 6.19(4) 


0.8(1) 


7(1) 1 0.0008 1 0.1 1 ~95% 1 


1 1.64 


0.1 1 6.8(3) 1 6.15(3) 


1.0(1) 


n.a. 1 0.0007 | 0.1 | ~ 65% | 



The first choice above (for n = l) was introduced by M. Liischer for conceptual purposes 
lITsll . and applied by Fukaya and Onogi in Schwinger model simulations l[l9Ll2flll . The 
question is whether one can conciliate the advantages mentioned above, with reasonable 
lattice sizes (say La ~ 1 — 2 fm), without increasing lattice artifacts and with a correct 
and reasonably decorrelated samplingof interesting observables. A first report of our 
ongoing study was presented in Ref. ll2lh . 

Results. Since gauge actions of the type (|1I2I3I| are non-linear in the link variables, 
the heat-bath algorithm cannot be applied. Instead we use a local HMC algorithm 
ll22] . which is competitive with heat-bath in the standard case. HMC trajectories have 
discretization dt quoted in Table[Tl and the trajectory length is 1. The volume is chosen 
to in order to allow a reliable determination of rg/a ~ 7 and in order have mild finite 
volume effects. We computed ro/a following a standard procedure lEsll . Our preliminary 
results are summarized in Table [Tl We estimated the topological charge with cooling and 
searching for the first plateau ll24ll . Since we cannot reliably measure the autocorrelation 
of the 2top, we quote - as an indicator of stability - the number of jumps of Qtop divided 
by the number of trajectories in the full history (//). Since we save a configuration 
only every 50 trajectories, the measured fj is only a lower bound on the frequency 
of jumps, which can be reliable only for fj <^ 0.01 (which is the interesting case for 
us). In particular for the Wilson action at /3 ~ 6.17 we expect fj ^ 0.015, which we 
measure. The stability of Qtop has to be compared with the autocorrelation of a typical 
observable (we quote the plaquette value under tP'^I). It is interesting to see that the latter 
is strongly decreased for non-zero £^ (at fixed ro/a). We have also studied actions of 
type (jSj), which have a smooth bound on the plaquette value, and therefore are better 
suited for efficient simulations with global HMC and dynamical fermions. Results will 
be presented elsewhere. We also checked that the results are consistent independently 
on the starting configuration. This is important because the constraint on the plaquette 
value could in principle generate more obstructions than the topological ones, and this 
would not be noticed simply from the autocorrelation. 

It has been pointed out [25] that the actions (H]) do not allow for the existence of 
a positive definite transfer matrix. However we checked that all the actions (I1I2I3I) 
have i'zYe-reflection positivity, which at least ensures the existence of a positive squared 
transfer matrix [26]. Moreover we have checked that the unphysical behaviour of the 
short distance force - which was observed in some cases, and related to the lack of a 



positive transfer matrix iEtIi - does not appear in our cases. 

We stress that - inspired by Ref. |Q] and by simplicity - our search is restricted to 
actions which are functions of the plaquette value, thus excluding other alternatives, 
which also affect the stability of Qtop, as for example in Ref. iEsLIiqIi . 

Conclusions. Topology conserving gauge actions could be highly profitable in QCD 
simulations. The suppression of small plaquette values may speed up the simulations 
with dynamical quarks. A stable 2top is useful in particular in the e-regime. We are 
investigating such actions, in view of the physical scale and the topological stability. 
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